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Abstract 

The heat quenching factor Q' (the ratio of the heat signals produced by nuclear and electron 
recoils of equal energy) of the heat-and-ionization germanium bolometers used by the EDEL- 
WEISS collaboration has been measured. It is explained how this factor affects the energy scale 
and the effective quenching factor observed in calibrations with neutron sources. This effective 
quenching effect is found to be equal to Q/Q' , where Q is the quenching factor of the ionization 
yield. To measure Q', a precise EDELWEISS measurement oi Q/Q' is combined with values of Q 
obtained from a review of all available measurements of this quantity in tagged neutron beam ex- 
periments. The systematic uncertainties associated with this method to evaluate Q' are discussed 
in detail. For recoil energies between 20 and 100 keV, the resulting heat quenching factor is Q' = 
0.91 lb 0.03 lb 0.04, where the two errors are the contributions from the Q and Q/Q' measurements, 
respectively. The present compilation of Q values and evaluation of Q' represent one of the most 
precise determinations of the absolute energy scale for any detector used in direct searches for dark 
matter. 

PACS numbers: 29.40.Wk, 95.35.-|-d 
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I. INTRODUCTION 



Heat-and-ionization germanium detectors are extensively used for direct search of the 
weakly interacting massive particles (WIMPs) that could constitute the Dark Matter halo 

MP on a Ge atom produces a nuclear recoil with 



of our Galaxy The scattering of a W 
a kinetic energy in the ten's of keV range 



i ^. The recoil is stopped in the detector volume 
within a distance of 100 nm. The energy loss occurs as a combination of ionization (elec- 
tronic dE/dx) and atomic collisions (nuclear dE/dx)|J|]. In heat-and-ionization detectors, 
the event is identified using two signatures. The first is the ionization signal, corresponding 
to the collection on electrodes of the electron-hole pairs created by the energy loss process. 
The second is the heat (or phonon) signal, recorded by a thermal sensor in contact with 
the germanium crystal. It has been demonstrated that combining the two simultaneous 
measurements provides an efficient discrimination against the large background of electron 
recoils originating from the natural 7 and (3 radioactivityP, Q]. The basis of this discrimina- 
tion is that the number of electron-hole pairs created by an electron recoil of a given energy 
is three to four times larger than that created by a nuclear recoil of the same energy. 

The process at the origin of the reduced ionization yield has been extensively studied since 
the original work of Lindhardj4|. The effect has been measured repeatedly as a function of 
recoil energy in germanium ionization detectors ^ liquid Nitrogen temperature (77 K) P, 

In these experiments, the ionization 



recoil energy in germanium ionization detectors at 
y, fl 3 y, Q, and more recently at 35 mK 



signals recorded with such detectors are first calibrated using gamma-ray sources, producing 
energetic electron recoils. This energy scale is called keV-equivalent-electron (keVge). Using 
this calibration, the detector is then exposed to a neutron source. The elastic collisions 
of neutrons with atoms in the detector volume produce nuclear recoils. The energy of 
the recoils is constrained by the use of monoenergetic neutron beams, the detection of the 
„eu.o. angle and/o. «.e-of-mg.. 1 fl fl 3 B Q • It . tKen observed tKat 

the ionization yield for a nuclear recoil is a factor Q ~ 0.25 smaller than that produced by 
an electron recoil of equal energy. 

In analogy with the ionization quenching factor Q, one can similarly introduce a quenching 
factor Q' for the heat signal in thermal detectors. This factor does not affect the ability 
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of heat-and-ionization detectors to discriminate electron and nuclear recoils^. However, 
this factor enters in the determination of the energy calibration for nuclear recoils and, 
consequently, in the energy threshold for their detection. Heat-and-ionization detectors 
are designed so that the initial deposited energy thermalizes in time for the signal read-out, 
independently of the process at its origin, with as little as possible losses to the outside world. 
It is thus expected that Q' is close to one. Possible deviations can only arise from processes 
that affect differently electron and nuclear recoils. In Ge crystal at cryogenic temperatures 
(typically 10 to 100 mK), the possible sources of such differences are usually considered 
to be small, but have never been measured precisely. Examples of such processes are the 
storage of energy in stable crystal defects generated by the recoiling nuclei, or the emission 
of photons during the initial stage of ionization. Systematic detector-dependent effects could 
also appear if the heat signal is read out before allowing for the full thermalization of the 
phonon excitations in the crystal. 

A direct measurement of Q' is a delicate experiment. Fig. Q summarizes the available 
measurements for all type of low-temperature devices. They are all consistent with unity, 
although the only measurement directly relevant to germanium, performed with a tagged 
neutron beam, is not very precise [11]. The measurements in the other bolometric devices 
involve ^oep]^ recoils from ^^''Pb alpha decays close to the surface of diamond and 
Te02 detectors. It could be argued that processes that could lead to deviations from 
unity in a germanium semiconductor might not show up as strongly, or even cancel out, 
in the other measured substrates. The technique of using surface recoils cannot be apphed 
easily to germanium heat-and-ionization detectors because of the presence of electrodes, 
dead layers and systematic effects in the drift of electrons and holes near the surface of the 
detector Q. 

However, as it will be shown here, heat-and-ionization bolometers provide an accurate 
measurement of the ratio Q /Q' . This idea has already been qualitatively stated in Ref. jlfil ]. 
where it was observed that the data appeared to be consistent with Q' = 1, without defining 
an explicit procedure for extracting a quantitative value. In this work, we define a procedure 
in which the Q' value for Ge recoils in Ge is obtained by combining the measurements of 



^ Except in the very special case where Q'=Q, clearly excluded by the measurements discussed in the 
following. 
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Q/Q' from heat-and-ionization detectors with the available direct Q measurements. In 
the first section, the mathematical formulae necessary to obtain the final result and, more 
importantly, its systematic uncertainty, are presented. In the second section, we review the 
available Q data, compute their average values as a function of recoil energy, and evaluate 
their statistical and systematic uncertainties. The third section presents the results oi Q/Q' 
measurements using the EDELWEISS detectors, and discusses their associated systematic 
errors. The results are also compared to those of Ref. In the fourth section, the Q/Q' 
measurements are compared with the computed average Q values. The final result can be 
interpreted either as the most precise available measurement of the heat quenching factor 
Q' for Ge atoms recoiling in the bulk of a Ge crystal, or, conversely, as an evaluation of the 
uncertainties on the absolute energy scale of nuclear recoils in the EDELWEISS heat-and- 
ionization detectors due to imperfection of the understanding of the underlying physics. 



II. DEFINITION OF TERMS 
A. Ionization measurement 



The amplitude of the ionization signal Ai is proportional to the number of collected 
electron-hole pairs A^. In heat-and-ionization detectors, the applied difference of potential 
on the electrodes V is chosen large enough so that losses due to trapping are small. This 
is verified^ by observing that Ai is independent of V . Under these conditions, for electron 
recoils 

E 

Aj^^ (X = — 

where E is the electron recoil energy and e-., is the average energy necessary to create 
one electron-hole pair. At 77 K, e^= 2.96 eVjlfi'|. At cryogenic temperature, this value is 
expected to increase slightly due to the variation of the gap energy. A value of 3.0 ±0.1 
eV is adopted here, consistent with the measurement of Ref. and the value used by the 
EDELWEISS and CDMS experiments. 

The amplitude Aj is calibrated using gamma-ray sources to provide Ej, the energy in 
units of keVee. 

2 In EDELWEISS detectors, this is verified for |y| > ~3 Volts. 
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For nuclear recoils, the average energy per electron-hole pair is approximately 12 eV, 
and varies with energy. The quenching factor for the ionization signal of nuclear recoils is 
defined as Q = e^/e„, and 

E E 

Aj^ri (xNn = — = Q — 

The Lindhard theory describing the energy losses of energetic ions in matter provides 
some predictions for the value of Q Within this model, the fraction of the incident 
energy dissipated in ionization is calculated from the electronic and nuclear stopping power 
dE/dx of Ge ions in Ge. Because of the approximations behind this calculation, it cannot 
be expected to provide an accurate prediction of the absolute value of Q. However, it 
will be shown in the following section that it describes reasonably well the overall energy 
dependence of Q. The Lindhard model will be used to interpolate between measurements of 
Q at different energies. We thus introduce this prediction as a reference function, Qref{E), 
evaluated for Ge recoils in Ge following the usual parametrizations P, 0] : 

where e is a dimensionless energy, k is related to the electronic energy loss [dE / dx electron = 
fc-^e), and for Ge recoils in Ge: 

e = ^E = 0.00354 E(keV) (2) 

^2/3 

k = 0.133— = 0.157 (3) 

^1/2 V / 

g(e) = 3e°-^^ + 0.7e°-^ + e (4) 



B. The heat or phonon measurement 

For electron recoils in heat-and-ionization detectors, the heat signal amplitude Ah is due 
to the initial energy of the recoil E and to the Joule heating associated to the current of 
collected electron and holes through the detector (Neganov-Luke effect jisl] ) : 

Ah,^ oc E + N^eV = {l + eV/e^)E (5) 

where e is the charge of the proton and V is the absolute value of the applied potential. 
For simplicity, we note the quantity eV/e^ as v, and An^-y oc (1 + v)E. For nuclear recoils in 
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a truly calorimetric detector {Q'=l), the amplitude is: 

AhAQ' = 1) « E + NneV = (1 + vQ)E 

If the phonon yield for nuclear recoils is reduced by a factor Q', the Joule heating contri- 
bution remains proportional to the number of created electron-hole pairs and: 

AH,n oc {Q' + vQ)E 

Once the heat signal amplitude An^-y is calibrated in units of keVee using a gamma-ray 
source, we obtain the energy Eh^^ = E and 

C. Recoil energy measurement for Q'=l and 

The two independent quantities E^ and Ej can be used to evaluate the following two 
quantities: 

X = {1 + v)Eh -vEi (6) 
y = Ej/x (7) 

We also define {x-y,y^) and (a;„,|/„) as the values of {x,y) in the case of electron and 
nuclear recoils, respectively. In the case where Q' = 1, the variable x and y have simple 
interpretations. As x-y = x„ = E, the variable x represents the recoil energy, irrespective 
of the type of recoil. The variable y measures the quenching effect relevant to the incident 
particle, as = 1 for electron recoils, while yn = Q for nuclear recoils. 

In the more general case where Q' 7^ 1, the relationships for electron recoils x^ = E and 
7/^ = 1 are still true. For nuclear recoils, however, the following values are obtained: 

Xn = Q'E 
Vn = Q/Q' 

A plot of the measured values of y versus x, as it is typically done|l[Q|, will still display 
the characteristic separation of the two populations of electron and nuclear recoils. However 
the position of the so-called "nuclear recoil band" {y{x) = yn{x)) has to be interpreted with 
care. First, the measured yn values will be divided by a factor Q' relative to the true value 



of Q. This has no consequence for the identification of nuclear recoils in Refs. , since 
the position of this band is always taken from the results of the neutron source calibrations 
of the relevant detector. However, the apparent recoil energy for nuclear recoil x„, obtained 
by assuming Q' = 1, is not equal to the true recoil energy E. For example, (5'=0.9, and an 
apparent recoil energy threshold of x„= 10 keV corresponds to a true threshold oi E = 11 
keV. 

The neutron calibration of the first heat-and-ionization detectors quickly established 
that Q' was not very different from 1, as the measured Un values were compatible to the 
direct measurement of Q performed at 77 K with Ge diode detectors. This work pursues 
this idea further by actually performing the evaluation of 

Q' = - (8) 

Vn 

using a compilation of the available direct measurements for Q and the most precise 
Un measurements performed by the EDELWEISS collaboration with its detectors. A com- 
plete assessment of the statistical and systematic uncertainties associated to both types of 
measurement will also be performed. 



D. Propagation of uncertainties 

To evaluate the systematic uncertainties on Q' as evaluated with Eq. |H1 it is useful to 
study its dependence on different variables. First, the uncertainty on Q' depends linearly 
on those on the Q measurement. In the next section, these uncertainties are discussed in 
details. As for the uncertainties on we have: 



dyn ^ {l+v)Q . s 

dEn Q'^E ^ ' 

dyn_ ^ {Q' + vQ) , . 

dEi Q'^E ^ ' 

dyn ^ QiQ'-Q) ( . 

dv g'2(l + v) ^ ^ 

An effect that must be taken into account when evaluating the ratio of the Q{E) and 
yn{x) measurements is that in general x ^ E. At lowest order (small deviations of Q' from 
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1 and small slope dQ/dE), the value of yn{E) can be obtained from that of yn{x) by: 

yniE)=yn{x) + ^^-^^E (12) 

where, as it will be shown later, dQ/dE can be taken from the Lindhard theory with a 
good precision. 



III. A REVIEW OF Q MEASUREMENTS 

Fig. 121 shows the results^ of the direct measurements of ionization quenching for germa- 
niun, .econ. in ge^aniu. fro. Refs. BflBay Q. TKe .eas.e.e.s .o,„ Ref. Q 

are discussed in Sect. I1V[ as they consist in a measurement of Q/Q' and not of Q. No large 
inconsistency can be observed. The recoil energy dependence of the data is well described 
by the Lindhard model (Eq. Q), shown as a continuous curve on the same plot. 

The measurements were performed with various experimental techniques. The systematic 
uncertainty on Q can be reduced by combining them. However, the recoil energy intervals 
used in the measurements are not the same, and one must take into account the energy 
dependence of Q when combining two measurements performed at different average recoil 
energies. For this reason, we will not perform weighted averages of Q{Er), but instead use 
Q{Eji)/Qref{Eji), with Qref{Eji) givcu by Eq. [TJ These weighted averages will be done for 
each of the recoil energy bin used in Sect. IIVI They are shown in Table lU In order to 
take into account possible correlations between the error bars within the same experimental 
data set in a given energy interval, only one measurement per experiment is used per inter- 
val. This measurement is chosen as the one with the smallest error. The weights are the 
square of the inverse of the quoted total uncertainties. The uncertainty on the combined 
Q{Eji)/Qref{Eji) ratios is computed assuming that the errors are not correlated between 
the measurements. If in a given energy bin, the reduced chi-square of the combination of 
N results, /{N — 1), is larger than one, the uncertainty is multiplied by the square root 
of this quantity. The average recoil energy < Er > for each energy bin is calculated as 
the energy of the individual measurements averaged using the same weights. Finally, the 
average ratio Q{Eji) /Qref{Eji) is multiplied with Qref{< >), yielding a set of values of 



^ The data from Ref. |a] are not included in this review because of the reasons explained in Ref. 
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Q{< En >) with their experimental uncertainties that exploits at best the available data 
sets. 

Table HI shows that the different measurements are consistent within small adjustment to 
the quoted errors, with four out of eight reduced x^values above unity, of which two are 
slightly above 2. This consistency is remarkable considering that the measurements were 
performed over a period of 38 years, with experimental setups and techniques that differ 
considerably and with constant improvements in the understanding of the systematic biases. 
Over the range from 5 to 200 keV, there is no sign of deviation of Q/Qref larger than ~±5%. 
This indicates that the choice of the function QrefiEn) for interpolating the results within 
an energy bin is appropriate. 

In the energy bin from 80 to 100 keV, the results of the three available measurements are 
consistent within less than 1%. The combined result has an uncertainty of 1%, making it one 
of the most precisely measured quenching factor in any type of detector. One can also note 
that the absolute value of Qref{ER) is very close to the experimental measurements. This 
is not too surprising since the overall normalizations of the nuclear and electronic stopping 
powers in the Lindhard model (Eqs. |21 to Ej) are based on experimental data. However the 
experimental values of Q{< E^ >) do not depend on this precise choice of reference. 



IV. THE EDELWEISS Q/Q' MEASUREMENT 

As discussed in Sect. Ill C\ the EDELWEISS heat-and-ionization detectors provide a mea- 
surement of y{x), where x = Q'E and, for events due to the elastic scattering of neutrons 
on Ge nuclei, the y{x) measurements have a Gaussian dispersion centered on the value 
Q/Q' . In this section, we present the results obtained with five detectors^, labeled GeAUO, 
GGAL GGA3, GSAl and GSA3. Their characteristics and performances are detailed in 
Ref. flQQ, where the e.pen,ne.ta> setup a.d ,ts ope.at.o. a.e also desenbed. The 
main difference between these detectors is whether the amorphous layer under the alu- 
minum electrodes is made of germanium (GGA series) or of silicon (GSA series), or whether 
there is none (GeAl series). This difference affects the charge collection properties of the 
detector P, The electrodes are polarized to a potential of 4.00 V, resulting in a value 



4 The two other EDELWEISS detectorsji^l GeA16 and GeA19 have not been used because of the stringent 
requirements for resolution and precision of the energy cahbration necessary for the present measurement. 
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ofv = 4/3. 

The experimental procedure and data analysis proceed according to the following steps. 
First, the detectors are calibrated in keVee using 7-ray sources (Sect. IIV A|) . They are 
then exposed to a fast neutron flux of the order of 1 to 10 n/cm^/h using a weak 252Qj? 
source. From this data is deduced i/n, the average value of ?/ as a function of x for elastic 
neutron scattering events fSect. HVB|) . At this point, the results are compared with those of 
Ref. {3]. The data are then corrected for multiple scattering and energy shifts fSect. HVT^ 
and the systematic uncertainties on the measurement are evaluated (Sect. IIVD|) . Finally, 
the quenching of the heat signal Q' is obtained by dividing the Q values of Sect. IIIII bv the 
corrected ijn values. 

A. Calibration with 7 rays 

The first step is the calibration of the detectors in keVee with 7-ray sources. An often 
overlooked difficulty in quenching factor measurements is that of the precise calibration of 
the detector in the low-energy range relevant for dark matter searches (10 to 30 keV, and 
even below for quenched signals). Eqs. 1^ and ITUl show that, for ~ 1 and f = 4/3, a 1% 
precision on ?/„ requires a precision of 1.3% and 0.7% on the measurement of the heat and 
ionization signals, respectively. 

It is difficult to calibrate the response of the bulk of a detector to ~10 keV 7 rays because 
of their short attenuation length, which excludes the possibility of using external sources. 
This difficulty is somehow alleviated in semi-conductor detectors. Because of the good 
linearity of their ionization signal, the extrapolation of the calibration from higher energy is 
more reliable than, for example, scintillating detectors. In the case of germanium detectors 
in an underground and very-low background environment, there is the additional unique 
possibility to check the calibration at 10 keVee using the doublet arising from the slow decay 
of the cosmogenic isotopes ^^Ge and ^^Zn (T1/2 = 271 and 244 days, respectively), and also 
the isotope ''^Ge arising from the activation of the detector with the 252 neutron source 
{Ti/2 = 11-5 days). Fig. El shows the two peaks at 8.98 and 10.34 keV coming from the 
electron conversion decay of Ge and Zn, respectively, as recorded in EDELWEISS detectors 
in a low-background run following a neutron calibration. With these data, it is possible to 
control with a precision of 1% the absolute energy calibration at 10 keVee of the ionization 
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and heat signals independently. At other energies, the relative calibration of the two signals 
can be easily checked using the distribution of the ratio of the two signals recorded in the 
Compton plateau produced by a ^'^^Cs 7-ray source. 



B. Neutron calibration data 

Fig. m shows the distribution of y{x) recorded in the detectors GSAl, GSA3 and GGA3 
when exposed to the ^^^Cf neutron source. The data of these three detectors are combined 
in that plot because they have similar experimental resolutions. 

Of the 20340 counts in Fig. ^ the overwhelming majority is part of four well-understood 
populations, described here. The first two populations are by far the most important. The 
first, centered a,t y = = 1, corresponds to 7 rays. It is well contained within the 90% 
efficiency zone computed from the experimental resolution on the heat and ionization chan- 
nels for these detectors, also shown on that figure. The second population lies mostly inside 
the zone expected for Ge nuclear recoils. On Fig. |5 the center of this zone is parametrized 
as in Ref. [2] as a band centered at y^ef = otx^ , with a = 0.16 and /3 = 0.18. Two lesser 
populations are observed, corresponding to the inelastic collision of a neutron on a ''^Ge 
nucleus with the excitation of the 13.3 and 68.8 keV states. In these cases, the value of y 
is the energy-weighted average of the quenching of the Ge recoil and of the electromagnetic 
radiation emitted in coincidence. Inelastic excitation of states decaying with higher-energy 
7 rays, observed as Compton electrons in the detector, would produce counts evenly spread 
between the two zones associated with electron and nuclear recoils, and over the entire 
energy range displayed in Fig. HJ 

In Fig. I3J there is a clear trend for y to increase with recoil energy for elastic nuclear 
scattering events. We define yn{,x) as the average y values for this population, as a function 
of X. The variations of yn{x) within the energy intervals defined in Table H] is not negligible 
at the level of precision required here. Accordingly, this energy dependence is removed by 
considering distributions of the variable defined as: 



where a and (3 are chosen such as to describe best the behavior of yn{x). Examples of 
D{x) distributions recorded in the neutron calibration runs for different intervals of recoil 





(13) 



1 — axf^ 
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energy x are shown in fig. El In practice, the standard values of a = 0.16 and 

13 = 0.18 are first used. The centroid Dn{x) of the D distributions for elastic scattering 
events is extracted for each bin of x. The Dn{x) values are transformed into values of ?/n(x) 
using the inverse of Eq. ^| The results for the sum of the five EDELWEISS detectors under 
study are listed in the third column of Table El They are also shown in Fig. (HI where for 
clarity is divided by the corresponding prediction of the Lindhard theory, Qref{x). 

The Unix) values are then least-square adjusted with the function a'x^' , yielding values of 
a' = 0.145 and f3' = 0.201. The entire procedure is then repeated, but this time using the 
new values of a' and /?' in Eq. El With these values, it is observed in Fig. that the new 
D{x) distributions are well centered on Dn = 0. However, the values extracted from 

the new distributions are equal to the previous ones, indicating that the procedure to remove 
the energy dependence of ?/„(a;) does not depend strongly on the precise choice of a and /3. 

In Fig. ini the EDELWEISS data are compared to those from Ref. ^|. For both data 
sets, only statistical errors are shown. Except for the values at x = 35 keV, the experimental 
results are compatible within the statistical uncertainties. However, two important sources 
of bias have not been corrected yet, and will be discussed in the following section. As Ref. 
was the first report of the observation of quenching effects in heat-and-ionization detector, it 
did not present the detailed information necessary for evaluating these biases. In particular, 
there is no quantitative statement concerning the precision of the 7-ray calibration, and no 
discussion of the effect of multiple scattering. As it will be shown in the following section, 
the multiple scattering of neutrons inside the detectors decreases the apparent value of ?/„. 
This effect could explain why the data from the ~180 g central volume of the EDELWEISS 
detectors seem to lie systematically below those obtained with the relatively more compact 
detector of Ref. hi 



C. Multiple scattering and energy scale corrections 

In order to interpret the measured values as Q{E)/Q'{E), the data must be cor- 

rected for two systematic effects: multiple scattering and the shift in energy scale when the 
variable x is transformed into a true recoil energy E. 

It was shown in Ref. [2^ that the fact that there are more than one neutron-nucleus 
interaction in a significant number of neutron scattering events has the consequence of 



13 



shifting down the measured values of y{x). For example, in an event where a single neutron 
produces two nuclear recoils of energy Ei and E2, the apparent quenching Qeff{E = E1 + E2) 
is equal to the energy- weighted sum of Q{Ei) and Q{E2). Depending on the relative values 
of El and E2, the weighted sum will lie somewhere between Q{E/2) and Q{E). As dQ/dE > 
0, this means that the measured value of Q{E) decreases as the contribution from multiple 
scattering increases. 

This bias can be corrected with the help of Monte Carlo simulations. In EDELWEISS, 



the results of simulations 



21[ based on the computer codes GEANT3 and MCNPX were 



compared and found to give consistent results to within 1%. The corresponding correction 
factors to be applied to are listed in Table ITTl They are of the order of 6 to 9%. They are 
larger by 2 to 3% compared with earlier simulations performed for Ref. [2^ with GEANT3. 
Since then, problems with these simulations have been identified and solved, leading to a 
better agreement with MCNPX. However, since this correction relies entirely on simulations, 
the uncertainty on the correction procedure is taken as a third of the correction. 

The second correction corresponds to the application of Eq. El It requires the knowledge 
of Q', and thus must be evaluated iteratively, first by evaluating it assuming Q'=l, and then 
replacing the derived Q' value in Eq. and repeat until convergence. Since Q' is very 
close to unity, only one iteration was found necessary to obtain stable corrections, listed in 
Table ITTl The size of the correction is approximately 2%. 

These two corrections are applied to the values of yn{x) in the third column of Table ITTl 
in order to obtain the Q{E)/Q'{E) values in the last column of the same table. 



D. Systematic uncertainties 

The seventh column of Table IHl lists the systematic uncertainties identified in addition 
to that associated with the multiple scattering correction. The sources of systematic biases 
considered here are those associated with the data selection and with the uncertainties on 
the energy calibration of Ej and E^, and on the value of v. 

The energy calibration is reliable to within 1%. In addition to the test described in 
Sect. IIV A| it was checked that the values of y^, the centroid of the y{x) distributions for 
electron recoils, stay within 1% of unity from run to run, whether with a neutron or 7-ray 
source or in the low-background runs. For the contribution of v, the dependence of the 
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heat gain as a function of the apphed voltage (Eq. was checked. Unaccounted deviations 
of more than 0.1 V are excluded, resulting in an uncertainty on ?/„ of less than 1%. The 
uncertainty on due to the uncertainty on e^(Sect. Ill A|) is less than 0.5%. 

Finally, the determination of the yn{x) values was repeated for each of the five detectors 
individually. This not only checks the effects of detector-to-detector fluctuations of the 
calibration, but also the influence of the different experimental cuts. For example, the lower 
cut on Ei is 2.5 keV for GGA3, GSAl and GSA3, and is 3.5 keV for GeAUO and GGAl. 
As di.cu.ed Ret. ^ correspond, to thresholds on recoil energy ot approximately 
11 and 14 keV, respectively. The effect of this cut can be observed on Figs. ID and El as a 
decrease of efficiency at the very lowest value of y{x). For this reason, for x values below 20 
keV, values of D{x) below -0.1 are excluded from the analysis. As a result, the fluctuations 
in yn{x) obtained from one detector to another for x < 20keV are larger than for most other 
intervals (Table [ij. 

For each interval in x, the largest deviation of any detector from the quoted yn{x) value 
is taken as a systematic uncertainty. This value is added in quadrature to 1% of the value 
of yn{x) to yield the systematic uncertainty listed in Table ITTl This conservative error tests 
possible deviations due to the technical differences between the GeAl, GGA and GSA detec- 
tors. It also tests the robustness of the measurement and in particular the reproducibility 
of the calibration procedure. 

This systematic error is added in quadrature with the uncertainty associated with the 
multiple scattering correction and the statistical error, resulting in the total error quoted in 
the last column of Table El These three contributions represent typically relative errors of 
the order of 3%, 3% and 0.5% of the quoted Q/Q' values. 

Two further systematic checks were performed. In the first one, the data selection and 
analysis procedure was applied to the simulated data sample described in Sect. IIV Gl The 
yn values extracted from the single scatter events are found to agree with the ax'^ parame- 
terization that had been input into the simulation, within a precision of 1%. The systematic 
biases due to the finite energy resolution and the selection procedure are thus well within 
the quoted total systematic errors. In the second test, values of yn were also extracted 
from inelastic scattering events. In this test, the data were corrected for the presence of an 
unquenched 7 ray with an energy of either 13.3 or 68.8 keV. The remaining ionization 
and heat energies. Eh — E^ and Ej — E^, were input to Eqs. (HI and [3 Within statistical 
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errors, the resulting y{x) distribution shows a band centered at the same values of yn{x) 
as the equivalent distribution for elastic neutron scattering events. Because of the reduced 
statistics (~1% of the elastic scattering sample), these data sets were not included in the 
final analysis. 

The present method to evaluate Q' assumes that the ionization quenching Q does not 
change between 17 mK (the temperatures at which Q/Q' is measured) and 77 K (the tem- 
perature for almost all of the Q measurements). The only direct experimental test of this 
assumption is provided by the Q measurements at 35 mK of Ref . ^| . They are compatible 
with those performed at 77 K, within an experimental uncertainty of approximatively 5%. 
However, there is no compelling reason to suggest that e^/e„ depends on temperature, since 
the few percent variation of the gap in germanium over this temperature range should affect 
equally both types of recoils. Therefore, the temperature dependence of Q between 17 mK 
and 77 K is assumed to be negligible. 



V. Q' RESULTS AND DISCUSSION 

Table UTTl summarizes the measurements of Q{E) and Q{E)/Q'{E) from Sect. Illll and IIV|, 
respectively. These two data sets are also compared in Fig. [7^. Because of the correction 
for multiple scattering and the relation between x and E, the best fit to the Q/Q' data 
(a = 0.149 and (3 = 0.209) differs slightly from the a' and (3' values describing the ?/„ data 
of Section FlV Bl This figure shows that the values of Q{E)/Q'{E) are systematically larger 
than those of Q{E). As a consequence, the resulting values of Q\E) deduced from their 
ratio are systematically lower than unity (see Table IIIII and Fig. [7|d). The weighted average 
of all measurements is 0.91. As the uncertainty on the Q' measurements at different energies 
are correlated, it is not possible to reduce the error on this average. However, it can be said 
that all values and errors between 20 and 100 keV are consistent with Q' = 0.91 ± 0.05. 
Outside this energy range, larger discrepancies are not excluded, as the measurements are 
less precise. Typically, the contributions of the Q{E) and Q{E)/Q'{E) measurements to the 
total uncertainty are ±0.03 and ±0.04, respectively. The latter error could be reduced by 
improving our understanding of the influence of multiple scattering. 

This determination of Q' is the most precise measurement of quenching in the bulk of 
a low-temperature bolometric detector (Refs. Q, plotted on Fig. [T]). It is three times 
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more precise than the direct germanium measurement )1 1[| . It is the first precise measurement 
based on data other than Pb and Po surface recoils from U/Th chain products. 

The measurement indicates that the thermal responses of our detectors to nuclear and 
electron recoils are different. This is an incitation to evaluate with more precision the effects 
that could be responsible. For example, simulations of the slowing down of 10 to 200 keV 
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germanium recoils in germanium with the simplest version of the program SRIM2003 
provide some indications that a few percent of the kinetic energy of the initial recoil can be 
trapped in defects due to displacements of atoms in the matrix. However the actual fraction 
depends considerably on the details of the simulation. How long these defects persist in the 
out-of-equilibrium environment along the recoil track is also a delicate question to address, 
especially in a low-temperature semiconductor medium 2^. 

It should be noted that the uncertainty related to the value of Q' does not affect the eval- 
uation of the efficiency for nuclear recoils, nor the rejection capabilities of the EDELWEISS 
detectors. This is because the bands for nuclear and electronic recoils are determined for 
every detector using a neutron source calibration, with exactly the same experimental condi- 
tions as in the physics runs. These bands are defined in terms of as observed in a given 
detector, and do not rely on any prediction concerning QiE) and Q'{E). However, the value 
of Q' does affect the interpretation of x in terms of absolute energy, and thus affects the 
overall efficiency of the detector because of the presence of an experimental threshold. The 
effect of a shift in energy scale on the limits on the WIMP-nucleon cross-section set by the 
EDELWEISS data was studied in Ref. jd]. In this case, a 10% shift of the energy scale was 
found to result in a 10% to 20% shift on the cross-section limits, a small effect at the scale 
of the present sensitivity of the experiments. In addition, the evaluation of both Q and Q' 
presented here are among the most precise quenching measurements for any detector used in 
the direct search for dark matter. In contrast, little is available concerning the scintillation 
quenching of iodine in Nal j^]. The experimental data on the quenching in Xe are only 
now just starting to converge and to cover the low energy range relevant for dark matter 
searches Q]. The wealth of available data on germanium and the consistency of the results 
make it one of the most reliable technique for detecting and identifying nuclear recoils. 
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VI. CONCLUSION 



The heat quenching factor (the ratio of the heat signals produced by nuclear and elec- 
tron recoils of equal energy) of the heat-and-ionization germanium bolometers used by the 
EDELWEISS collaboration has been measured. It is shown that the calibration of these 
detectors with neutrons provides a measurement of the ratio Q/Q', where Q and Q' are the 
quenching factors of the ionization yield and heat measurement, respectively. Consequently, 
the existing direct measurements of Q have been reviewed and the resulting values have been 
divided by the EDELWEISS Q/Q' measurements. The resulting heat quenching factor for 
germanium recoil energies between 20 and 100 keV is Q' — 0.91 ± 0.03 ± 0.04, where the two 
errors are the contributions from the Q and Q/Q' measurements, respectively. The evalua- 
tion of both Q and Q' presented here are among the most precise quenching measurements 
for any detector used in the direct search for dark matter. 
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(keV) 


Refs. 


Q / Qref 




<Er> 

(keV) 


Qrefi< Er >) 


Qi< Er >) 


5-10 


\9] 


1.05 ± 0.08 




8.6 


0.230 


0.241 ± 0.019 


10-15 


\7, 9] 


0.99 ± 0.10 


2.01 


12.2 


0.241 


0.239 ± 0.024 


15-20 


^7, 9] 


0.95 ± 0.04 


0.34 


16.4 


0.252 


0.240 ± 0.011 


20-30 




0.98 ± 0.03 


1.29 


24.9 


0.269 


0.265 ± 0.009 


30-40 


[5, 9] 


0.98 ± 0.03 


1.34 


35.4 


0.284 


0.279 ± 0.008 


40-60 




1.05 ± 0.03 


0.48 


54.9 


0.306 


0.322 ± 0.009 


60-80 




1.03 ± 0.02 


0.71 


74.0 


0.323 


0.331 ± 0.007 


80-100 




1.02 ± 0.01 


0.31 


89.3 


0.334 


0.342 ± 0.005 


100-150 


[5, 7,iO,i]J 


1.00 ± 0.02 


2.21 


113.6 


0.350 


0.349 ± 0.008 


150-200 


\W] 


1.01 ± 0.02 




178.2 


0.384 


0.388 ± 0.007 



Table I: Average of measurements of Q, the ionization quenching of Ge recoils in Ge, as a function 
of the recoil energy Er. 



Er 


<Er> 


yn{x) 


Stat. 


Multiple scattering 


E 


Calib. 


Q/Q' 


(keV) 


(keV) 




error 


correction 


(Eq. 12) 


error 




10-15 


12.8 


0.226 


0.014 


1.039 ± 0.013 


0.985 


0.017 


0.231 ± 0.077 


15-20 


17.4 


0.252 


0.003 


1.059 ± 0.020 


0.983 


0.008 


0.263 ± 0.029 


20-30 


24.4 


0.276 


0.001 


1.070 ± 0.023 


0.983 


0.003 


0.290 ± 0.012 


30-40 


34.4 


0.295 


0.001 


1.083 ± 0.028 


0.982 


0.004 


0.314 ± 0.014 


40-60 


48.4 


0.315 


0.001 


1.096 ± 0.032 


0.982 


0.003 


0.339 ± 0.013 


60-80 


68.7 


0.341 


0.002 


1.088 ± 0.029 


0.981 


0.005 


0.364 ± 0.017 


80-100 


89.1 


0.357 


0.003 


1.084 ± 0.028 


0.980 


0.004 


0.380 ± 0.014 


100-150 


119.3 


0.375 


0.003 


1.083 ± 0.028 


0.979 


0.006 


0.398 ± 0.016 


150-200 


170.1 


0.391 


0.006 


1.057 ± 0.019 


0.978 


0.012 


0.404 ± 0.031 



Table II: Measurements of Q/Q' for germanium recoils in germanium obtained from the neutron 
cahbration of the EDELWEISS detectors GeAUO, GGAl, GGA3, GSAl and GSA3. 



20 



(keV) 


<Er> 

(keV) 


Q{E)/Q'{E) 


Q 


Q' 


10-15 


12.8 


0.231 ± 0.077 


0.239 ± 0.024 


1.04 ± 0.36 


15-20 


17.1 


0.263 ± 0.029 


0.2-10 ± 0.011 


0.91 ± 0.11 


20-30 


24.4 


0.290 ± 0.012 


0.265 ± 0.009 


0.91 ± 0.05 


30-40 


34,4 


0.314 ± 0.014 


0.279 ± 0.008 


0.89 ± 0.05 


40-60 


48.4 


0.339 ± 0.013 


0.322 ± 0.009 


0.95 ± 0.05 


60-80 


68.7 


0.364 ± 0.017 


0.331 ± 0.007 


0.91 ± 0.05 


80-100 


89.1 


0.380 ± 0.014 


0.342 ± 0.005 


0.90 ± 0.03 


100-150 


119.3 


0.398 ± 0.016 


0.349 ± 0.008 


0.88 ± 0.04 


150-200 


170.1 


0.404 ± 0.031 


0.388 ± 0.007 


0.96 ± 0.07 



Table III: Summary of the results of the compilation of Q measurements, of the EDELWEISS Q/Q' 
measurement, and the deduced values of Q' as a function of recoil energy. 
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Zhou 94 



Rb on C 



Alessandrello 96 



Alessandrello 97 



Simon 03 



Pb on Te02 



Pb/Po ot\ Te02 



Ge in Ge 



0.8 



0.9 



Figure 1: Available measurements of Q', the quenching of the heat sig nal for nuclear recoils, in 
germanium, diamond and Te02- The measurements are from Refs. 12 1 (Zhou 94), (Alessan- 
drello 96 and 97) and [ll| 
100 keV. 



(Simon 03). In all cases, the recoil kinetic energies are approximately 
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Figure 2: Experimental results of the direct measurement of the ionization quenching for germa- 
nium recoils in germanium, from Refs. ^ (Chasman 65), 3| (Chasman 67 and 68), a] (Jones 71 
and 75), ,9] (Messous 95), (Baudis 98) and (Simon 03). The line represents Eq. Q with 
parameter values as of Eqs.|21to|2 
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EDELWEISS calibration at 10 keV 
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Figure 3: Example experimental checks of the calibration at 10 keVge of the EDELWEISS detectors 
GSAl (top) and GGA3 (bottom). The full lines indicate the expected position of the activation 
peaks at 8.98 and 10.34 keV. Left: heat signal spectra. The FWHM resolution of the peaks are 
0.76ib0.05 and 0.58ib0.05 keV for GSAl and GGA3, respectively. Right: Ionization spectra. To 
improve the identification of the two peaks, these spectra were obtained using gates on the 8.98 
and 10.34 peaks in the heat spectra on the left. The FWHM resolution of the ionization peaks are 
1.3ib0.1 keV for both detectors. 
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EDELWEISS ^^Xf Calibration 




X = Recoil Energy (keV) x Q' 



Figure 4: Experimental data recorded with a Cf neutron source with the EDELWEISS detectors 
GSA3, GSAl and GGA3. 
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Figure 5: Experimental distributions of D = {y — ax^)/{l — ax^) for different values of recoil 
energy interval x, for the sum of the five EDELWEISS detectors used in the present work. Here, 
a = 0.145 and (3 = 0.201. 
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Figure 6: Ratio of the experimental values of y„ divided by the Lindhard model for EDELWEISS 
(full circles) and the data from Ref. The same ratio is also shown for the parameterization 

l/n = 0.145x''''^*'^, corresponding to the best fit to the EDELWEISS y„ data (dashed line). 
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Figure 7: (a): values of Q/Q' measured by EDELWEISS (full circles) and values of Q obtained in 
the compilation described in Sect. ITTll (open circles). The horizontal error bars correspond to the 
range of the different energy intervals used in the Q and Q/Q' analyses. The dot-dashed curve is 
the Lindhard parameterization described in the text. The dotted curve represents the best fit to 
the EDELWEISS Q/Q' data {Q/Q' = O.U9{ERf-'^°'^). (b): Q' values obtained from the ratio of 
these data. The dashed line is the weighted average of all values {Q' = 0.91). 
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